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We introduce the Lagrangian for a multi-scalar field configuration in a A''-dimensional internal 
space endowed with a constant metric Qi^. and generalize the quintom cosmological scenario. We 
find the energy momentum tensor of the model and show that the set of dual transformations, that 
preserve the form of the Einstein equations in the Friedmann-Robertson- Walker (FRW) cosmology, 
is enlarged. We show that the stability of the power law solutions leads to an exponential potential 
which is invariant under linear transformations in the internal space. Moreover, we obtain the general 
exact solution of the Einstein-Klein-Gordon equations for that potential. There exist solutions that 
cross the phantom divide and solutions that blow up at a finite time, exhibiting a superaccelerated 
behavior and ending in a big rip. We show that the quintom model with a separable potential 
can be identified with a mixture of several fiuids. This framework includes the ACDM model, a 
bouncing model, and a setting sourced by a cosmic string network plus a cosmological constant. 
The we concentrate on the case where the dimension of the internal quintessence sector A'^, exceeds 
the dimension of the internal phantom sector Nph- For {Nq,Nph) ~ (2, 1) the dark energy density 
derived from the 3-scalar field crosses the phantom divide and its negative component plays the role 
of the negative part of a classical Dirac Field. 
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I. INTRODUCTION 



A large number of recent cosmological data, including 
Type la supernovae [H , Large Scale Structure 0, [1] and 
Cosmic Microwave Background anisotropies [3| have pro- 
vided strong evidences for a spatially flat and accelerated 
expanding universe at the present time. For general rela- 
tivistic FRW cosmologies, such an accelerated expansion 
implies the existence of a mysterious exotic matter with 
large enough negative pressure, the dark energy, whose 
energy density has been dominating the recent stage of 
the universe. The astrophysical feature of dark energy 
is that it remains unclustered at all scales where gravi- 
tational clustering of baryons and nonbaryonic cold dark 
matter can be seen. The combined analysis of the dif- 
ferent cosmological observations suggests that the uni- 
verse consists of about 70% dark energy, 30% dust mat- 
ter (composed by cold dark matter and baryons), and 
negligible radiation. Although the nature and the ori- 
gin of dark energy are unknown, several candidates have 
been proposed to describe it. The simplest theoretical 
candidate is the cosmological constant A [H, @ whose 
equation of state is p = —p. However, this proposal suf- 
fers from the well-known cosmological problem according 
to which there is a unexplained extraordinary discrep- 
ancy (of about 120 orders of magnitude) between the 
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observed value of A and the value predicted by quantum 
field theory Q, [I]- For this reason, alternatives routes 
have been suggested to study the dark energy and its 
dynamical evolution So far a wide variety of scalar- 
field dark energy models have been studied, including 
quintessence Q, /c-essence [T^], tachyon [ll|, phantom 
[l^ . ghost condensate [ll] and more recently, quintom 
dark energy [13, [13, HI] • In addition, other possibihtics 
include interacting dark energy [Tst . braneworld models 

16 1 , Chaplygin gas models [T3| , holographic dark energy 

18| and many others. 

The quintom model is a hybrid with a quintessence 
component, usually constructed by a scalar field and a 
phantom scalar field minimally coupled to gravity. This 
model has aroused interest due to a debate about the 
possibility that the dark energy equation of state crossed 
the phantom divide in the near past. Notice that in the 
quintessence models the kinetic energy is always positive 
definite for a positive definite potential, whereas one has 
exactly the opposite in phantom models. So, neither the 
quintessence nor the phantom alone can fulfill the tran- 
sition a non-phantom to phantom stage. 

Many authors have studied the quintom paradigm with 
several potentials from the dynamical autonomous sys- 
tem perspective [m,[13- In this context, some of them 
have performed a general analysis by including the con- 
tribution of both a quintom component (dark energy) 
and a matter component into the Friedmann equation. 
For example, in Ref.p^ it was found that the phantom 
attractor is not generic, which means that in a quintom 
model there may exist either de Sitter attractors asso- 
ciated with the saddle points of the potential, or track- 
ing attractors in the asymptotic regime where the scalar 
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fields diverge. 

Finally, we want to stress that motivated by string the- 
ory a quintom scenario has been found by taking into ac- 
count the non-perturbative effects of a generalized Dirac- 
Born-Infeld action [g^l- Moreover, a "No-Go Theorem" 
has been proven to constrain the building of quintom sce- 
nario [i^l ■ According to this theorem in order to achieve 
the crossing of the phantom divide it is necessary to have 
more than one degree of freedom. Clearly, this theorem 
offers a concrete theoretical justification for the quintom 
paradigm. 

In this paper we present a multi-scalar field model em- 
bedded in a A'^-dimensional internal space structure en- 
dowed with a constant internal metric Qik- In Sect. II, 
we introduce the Lagrangian for the multi-scalar field 
model, calculate the energy momentum tensor and ob- 
tain the energy density and pressure of the multi-scalar 
field. We also study how the dual transformation, that 
links contracting and expanding cosmologies is realized 
when the internal structure is taken into account, and 
we then perform the stability analysis for power law so- 
lutions. In Sect. Ill , we report the general exact solution 
of the Einstein-Klein-Gordon equations for the exponen- 
tial potential and illustrate the existence of the crossing 
of the phantom divide, showing that the model has a 
smooth transition from the quintessence to the phantom 
era. In Sect. IV we study quintom and a 3-dimcnsional 
scalar field models driven by separable potentials. Fi- 
nally, in Sec. V the conclusions are stated. Throughout 
this paper we shall use natural units (G = c = 1). 



II. MULTI-SCALAR FIELD COSMOLOGIES 

We wish to construct a general class of dark energy 
models with late-time phantom behavior motivated by 
the fact that this peculiar sort of expansion fits confort- 
ably within current observational results. Accelerated 
expansion cannot extend indefinitely back in time, be- 
cause it would spoil the formation of structures, and well 
founded theoretical arguments support the view that our 
universe was dark matter dominated for a long time span 
before it entered the current dark energy dominated ac- 
celerated phase. In addition, results point out in the di- 
rection that dark energy at present does not just produce 
acceleration but what is rather called super-acceleration, 
thus indicating that dark energy is phantom-like, and it 
is quite likely that it was it actually became phantom- 
like just recently, i.e. dark energy may have crossed the 
so-called phantom divide. 

The mainstream approach (for simplicity and because 
there is no evidence in favor of the contrary) is to consider 
that dark energy and dark matter interact only gravita- 
tionally, so they "see" each other by their effects in the 
expansion. Intuition suggests, then, that the ability for a 
dark energy fiuid to cross the divide is not ruined by the 
presence of dark matter. For that reason, even if one's 
goal is to try and construct cosmological scenarios con- 



taining dark matter and dark energy featuring a crossing 
of the phantom divide, it makes sense to start off by con- 
sidering the dark energy fluid alone. Put another way, 
it is very likely that the construction of a dark energy 
source with the desired property leads when combined 
with dark matter to a cosmological scenario in which the 
dark energy bit retains the same feature. 

Quite a few scalar field model configurations with the 
ability of crossing the phantom divide can be found in the 
literature, and perhaps the most populated class is that 
of the generically dubbed quintom models [3, [13, HI]. 
Their main feature is that the kinetic energy is not a defi- 
nite positive quantity. Quintom scenarios have attracted 
considerable attention because the mentioned property 
of the kinetic energy allows giving a unified treatment of 
phantom and conventional scalar fields. 

Taking quintom models as an inspiration we assume a 
A'^-dimensional internal space structure endowed with a 
constant metric Qik and propose a new theoretical setup 
with multiple fields where i = 1...N. This compos- 

ite field can be thought of as a vector within the internal 
space. For any two vectors Ui and Wk belonging to this 
space their scalar product is u • w = QikUiWk and we 
denote the norm as w • w = = QikWiWk- In what 
follows we use latin indexes for the internal space and 
summation over repeated indexes will be assumed. 

The Lagrangian for this model is given by 

C = ^Qtkdt,Mx)d''(l)k{x) - V{(f>i,<j>2,-,(I)N), (1) 

where 1^(01,02, ■■■t4>n) is a selfinteracting scalar poten- 
tial depending on the fields (j)i,4'2, ■■■i4>n only. The en- 
ergy momentum T^j^ associated to the multi-scalar field 
Lagrangian ^ is obtained through the standard equa- 
tion = 25C/5g^^'' — g^uC^ so we get 



:Q^k^^cj),^^(Pk - V 



(2) 



Our framework is that of spatially fiat FRW cosmolo- 
gies and the energy density and pressure of the dark en- 
ergy fiuid are 

P = ^Qik'Pi^k + V, ]^Qik4>i4>k ~ V, (3) 

where we have assumed homogeneous scalar fields (f>i ~ 
4>i{t). The evolution of a universe sourced by that config- 
uration of multiple interacting fields in the fashion above 
is governed by the equations 



3i?' = ^Q^k'i)^'i}k + V, 



and they get combined to give 
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H = --Qik(t>i4>k- 



(4) 
(5) 

(6) 
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Our proposal is to use this extension which, due to the 
presence of an internal space, allows to consider a richer 
structure than in the precursor quintom cosmologies. 
The usual quintom model fits in our description for the 
choice Qik = diag (1,-1). In particular the internal 
metrics Qik = diag (1,0) and Qik = diag (0,-1) allow 
to reproduce to quintessence and phantom cosmologies 
respectively. 



A. Duality 

Interestingly, the internal space structure also allows 
to interpret duality transformations linking contracting 
and expanding conventional cosmologies in a very intu- 
itive and elegant fashion. In fact, the a ^ 1/a duality 
transformation implies H —H, and from the Einstein 
equations ^ and (O it follows that there are two al- 
ternative ways to achieve the transformation that leaves 
those equations unaltered: 



• Case 1 



(7) 



• Case 2 

Q^k ^ ~Q,k, V ^ Q,k^^<pk + V{(f>,) (8) 

A realization of the duality as given by case 1) requires 
potentials which remain real under the transformation of 
the fields (j)k ii'k- However, the realization induced 
by the internal structure and given by case 2) is of wider 
applicability as the former requirement on the potential 
is not necessary. In addition, both the transformed and 
non-transformed fields arc real. For instance, phantom 
cosmologies characterized by a future singularity occur- 
ring in finite time can be obtained by performing the 
a ^ 1/a duality transformation onto cosmologies with a 
final big crunch. For the case 1), one of them (phantom 
or big crunch cosmology) will be associated with imag- 
inary fields. However, for the case 2) both cosmologies 
are described by real fields. 

B. Stability of power law solutions 

We are going to investigate now the existence of 
asymptotic power law solutions by means of a structural 
stability analysis, in other words, we are going to unveil 
the asymptotic nature of potentials allowing such stable 
solutions. To this end we introduce the barotropic index 
7 = (p + p)/p, and from the definitions above and the 
Einstein equations (jH) and ^ one gets the dynamical 
equations for 7 = —2H/3H^, namely 



7 - (7 - 2) 3i77 + 



V 
V 



(9) 



Clearly, the solution 7 = 2 represents an equilibrium 
point, but it corresponds to non-accelerated expansion, 
specifically, a cx t^^"^. Given that practically any other 
interesting cosmological model one can think of will have 
7 < 2, let us demand then that Eq. ([9]) admits an- 
other equilibrium point with 7 = 70 representing in con- 
sequence a solution with a cx t^/^^". In particular, a so- 
lution representing a phantom late-time attractor would 
be characterized by a constant and negative value of 7. 
Below, we will show that the requirement of the existence 
of such an equilibrium point restricts the functional form 
of the V{(l)i) potential. Then we formulate a structural 
stability analysis of the equation governing the evolution 
of the parameter 7. This is equivalent to imposing the 
asymptotic condition on the potential V + SjoHV = 0, 
which can be integrated to give V = Voa~^'^° with Vq a 
positive integration constant. By choosing this potential 
the asymptotic regime of 7 is governed by the equation 



7 = 3i/(7-2)(7-7o). 



(10) 



Let us now define 7 = 2 — e and recast Eq. (fTO|) in 
terms of this new definition. By expanding the rhs of this 
equation up to order e, for expanding universes {H > 0), 
the solution 7 = 2 is a repcller. A similar analysis can 
be done for the case 7 = 70 7^ (the particular case 
70 = will investigated at the end of this section). If one 
defines 7 = 70 — e it turns out that the 7 = 70 solution is 
an attractor provided 70 < 2, which is on the other hand 
our working hypothesis. If one then calculates V and H 
as functions of t for a = i^/STo (^an be inferred from 
Eqs. (HI) and © that 



t2 ' 



2(2 -70) 



(11) 



As it has been shown, in spite of having started from 
a general formulation, we have been able to calculate the 
value of Vq and the asymptotic dependence of the poten- 
tial on cosmic time. However, in what the reconstruction 
of the potential is concerned we need to make contact 
with standard formulations in physical cosmology and to 
specify the kinetic energy as a function of the fields and 
their derivatives. 

For the attractor solution (which satisfies a power law) 
it can be seen that a satisfactory choice for the asymp- 
totic behavior of the scalar fields is 4> = (pQ log t with cf>Q 
a constant vector in the internal space. Besides, we are 
considering that physical quantities such as the potential 
energy depend only on scalar quantities in this internal 
space, so for simplicity we will consider the potential de- 
pends on the scalar fields through a linear form of the 
fields. To that end we introduce a constant vector a in 
the internal space. Under these assumptions we build the 
scalar quantity <p = cy. ■ (f>Q logt, which leads to 



V = Foe" 



-a (f) 



(12) 



where we have fixed a ■ (pQ — 2. One can also find the 
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following relations 



2a 



(13) 



Finally, in terms of the parameters of our framework the 
attractor solution gets recast as 



2/c 



(14) 



Due to the internal space structure of our model the ex- 
ponent of the latter can take either sign, so both phantom 
and non-phantom power law solutions are available. 

Finally, we conclude this section with the investigation 
of the stability of the solution 70 = 0. If we impose the 
fulfillment of Eq. ^ asymptotically one obtains V/V = 
with V = Vq a constant. Now defining 7 = — e and 
using again PTT)) it follows that de Sitter solution {H ~ 
Hq) is asymptotically stable. 



III. THE EXPONENTIAL POTENTIAL 

We have investigated the stability of the power law so- 
lution with a exponential potential and the possibility of 
phantom and no phantom behavior for the scale factor. 
In this section we shall find the general exact solution for 
our A'^-scalar field model driven by the exponential poten- 
tial in the spatially flat FRW spacetime. In this case 
the Eqs. (H])-® represent a non- linear coupled system of 
(N + l) differential equations for the scale factor a{t) and 
the N scalar fields 4'i{t). For a non-phantom scalar field 
driven by a exponential potential and a free scalar field, 
the general solution of the system of equations (H))-® 
was found in Ref. [l^. Using the Einstein equation ^ 
it is easy to prove that the geometrical object 



bk = akH + 

a"* 



(15) 



with Cfc a constant vector in the internal space, is a first 
integral of Eq. ([S]) . Inserting the latter into the Eq. ^ 
and expanding, we obtain 



-2H ^ a^H^ + 2a c 
Making the change of variables 



H c2 



-3/n 



T = a ■ c n 



(16) 



(17) 



in Eq. p6p . it becomes a nonlinear second order differ- 
ential equation for s(t): 



1 



4 cos^ 9 



5^"+! = 0, (18) 



where the dot stands for differentiation with respect to 
the argument of the function and 



cos( 



ac 



(19) 



The sign of the parameter n depends exclusively of the in- 
ternal metric Qik- In particular, for an Euclidean metric 
this parameter becomes negative. After the new variable 
s is calculated, we obtain the scale factor a in terms of 
the cosmic time t by using Eqs. pT]) . 

Inserting the first integral ([TS]) into the Friedmann 
equation ([4]), we get a quadratic equation in the expan- 
sion rate H which has real solutions only when its dis- 
criminant is positive definite, this condition can be recast 
as (a • c)^ -|- (6 — a^){c'^ + 2a^V) > 0. It gives a relation 
among the integration constants and the exponent of the 
exponential potential. 

Eq. (jlSp can be linearized and solved using the form 
invariance group of this nonlinear differential equation 
under a non local change of variable. This technique re- 
lates a class of nonlinear differential equations with the 
damped harmonic oscillator one. In fact, Eq. be- 
longs to the following general class of second order non- 
linear ordinary differential equations: 

s + n F{s)s + V F{s) [ F{s)ds + S F{s) = 0, (20) 



where s = s{t), F{s) is a real arbitrary function, and /z, 
i/, S are constant parameters. 

By using the new pair of variables {h^if}, defined by 
the non-local transformation 

F{s)ds = G{h)dh, F{s)dt = G{h)dr], (21) 

Eq. (|20p becomes the nonlinear ordinary differential equa- 
tion 

h" + ^iG{h)h' + vG{h) [ G{h)dh + 5 G{h) = 0, (22) 



where ' denotes differentiation with respect to 77. Inter- 
estingly, Eqs. ((20|) and ([22|) arc related between them by 
means of the formal change F ^ G and s h. Then, the 
non-local transformation pT|) preserves the form of these 
equations mapping solutions of Eq. (|20p into solutions 
of Eq. (|22p for any function F{s) and G(7/), thus linking 
solutions of two different physical configurations. Taking 
into account the nature of this non-local change of vari- 
ables it is not always possible to find explicit solutions. 
However, we can use the form invariance to relate solu- 
tions of a nonlinear equation (|20p with a linear one. In 
our case, choosing the specific parametrization F{s) = s" 
and G{h) = 1 in the change of variables (PT|) 



s"(is, 77 



(23) 



and introducing these new variables in Eqs. (|20p and 
I, they become 



^ n+l 

h" + iih' + vh + 5 ^0. 



(24) 
(25) 
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It shows that the non-local change of variable ((23)) trans- 
forms the non-linear differential equation into the 
linear damped harmonic oscillator equation ([25]) . 

Comparing our equation with Eqs. and (PS)) . 
for = 1, we have two different cases to analyze accord- 
ing to the n value, 

• n^-l and ^ = ^ F{s) = s\ 

• n = -1 and I. = ^ = i, 5 = 

Below we find the general solution of the Einstein-Klein- 
Gordon equations for the exponential potential. 



A. Explicit solution for n 7^ 1, —2 and v 



In this case explicit r-dependent solutions can be found 
introducing the function v[t) (for more detail see Ref. 
[2^) and making the substitution 



(26) 



n J ki + J v^^dr ' 



into the equation p^ . This reduces it to i) = 0, whose 
solution is v{t) = k2 + k^r. After inserting the solution v 
in Eq. ipBl) . integrating and using s{a) as given by Eg. lfTT)) . 
we obtain the scale factor 



aiv 



(n-f- l)(n + 2) 



1/3 



k\v\ 



-«ii/3 



(27) 



where the constant k has been expressed in terms of the 
old one [2^. Finally, we can get the components 0' of the 
multi-scalar field 4> integrating the Eq. (jisp 



dv 



dlna c^a 



dv 



a ■ c 



(28) 



To this end, we rewrite the second term of the Eq. ([28 
using the definition given in Eq. (jl7p as 



2 d 



dv 



In ki 



v''dT 



(29) 



Inserting the latter into the Eq. ([28]) we obtain the {v) 



:ln| 



jn + 2)c, 
n{a ■ c) 



In 



(30) 



where A(j) i(v) = (f>i{v) — (f>Qi and (j)Qi are N integration 
constants [30|. For positive n , when v —> 0, the leading 
terms in the scale factor is given by 



fcl' 



(31) 



So, the scale factor blows up near the origin v = and 
the components of the multi-scalar field have a logarith- 
mic divergence A0j cx In \v\ for all n. We find a resem- 
blance with the big rip behavior in the phantom universe 



reported recently by Wei |26j with an inverse power po- 
tential V = Vo(j)^^. 

As we have integrated the second order differential for 
the scale factor instead of the original Friedmann 
equation ([4]), we have to constrain the integration con- 
stants. By using Eqs.dH) and we obtain 3H'^+H = V, 
which can be rewritten in terms of the variable u = as 



d u 
dv^ 



3uV 
(a • c)2 



(32) 



After replacing the solutions u{v) and 4>{v) in Eq.([32|) the 
dependence in the variable v is missed and we obtain the 
following relation between the integration constants 



k = 



2c2|0o| cos^ 



sign 



(33) 



Hence, the scale factor p7p and the multi-scalar field ([30]) 
whose integration constants verify the relation p3p are 
exact solutions of the Einstein-Klein-Gordon equations 
O-®. 



B. Explicit solution for n = 1 and v = Q 

In the n = —1 case, following the same steps we get 
a=[\v\\h^+\n\v\\f\ (34) 

= -i In |5i + In \v\\ - — ^ In |&i + In (35) 

3 OL ■ c 

where bi is an integration constant. Inserting Eqs. 

into the Eq. ^ we find the relation = 

2|(/)o|c^ cos^ 6*, where = 6 and the constant has 
been renamed properly. 



C. Implicit solution 

The general implicit solution of Eq. ([TS]) can be found 
for arbitrary choices of the parameters by solving Eq. 
((25)) and using Eqs. (fT7|) . ((23)) . In this case we have two 
different set of solutions according to the value of n. For 
n —1, (5 = and yU = 1 the scale factor is given by 

-n/3(n+l) 



(77)= y/{n+ 1) {bi exp X^T] + &2 cxp A+77) 

(36) 

where A± are the roots of the characteristic polynomial 
of the Eq. ([25]), A± = (-1 ± ^/l-4l^)/2. Also, from Eq. 
rst we calculate the multi-scalar field 



a • c 



-V, 



(37) 



where (f>iQ are integration constants. Now, inserting Eqs. 
((5S)) and ((57)) into the Friedmann equation ((5^ the inte- 
gration constants are constrained to 



6162 



Vo 



2(1 -4i/)c2|0o|cos2 6'' 



(38) 
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Finally, in the n = —1, ly ~ and S = (4cos^ 9) ^ case 
we have 



a{r]) = exp 



a ■ c 



bi 



2|0o|c2 cos2 6l' 



(39) 
(40) 
(41) 



For 6 = Owe reobtain the standard formula for two scalar 
fields driven by an exponential potential Ref. [1^. 



IV. 2, 3-DIMENSIONAL INTERNAL SPACE 

Below we investigate several interesting examples with 
different internal metrics. 



A. Big Rip solution 

Let us explore a new class of quintom model associated 
with the internal metric 



iik — 



1 

1 



(42) 



As the sign of the kinetic term is non definite, a new 
feature arises, the possibility of having a smooth tran- 
sition between quintessence and phantom scenarios (see 
Fig.l). For < the scale factor has two branches. In 



6 

a(T)4 
2 




-2-10 1 2 3 

AT 

Figure 1: The two branches of the scale factor for — —1. 

the branch At < the universe expands and has a finite 
time span, it begins to evolve from an initial singularity, 
has a transient phase to a superacelerated expanding sce- 
nario ending in a big rip. Initially the universe is dom- 
inated by an quintessence scenario with kinetic energy 
T{4>ii4>2) > and barotropic index j = {p + p)/p > 0, 
after that, the universe evolves to a dominated phantom 
dark energy stage with T{(f)i, ^2) < and 7 < until the 



big rip occurs when At 0~ and a 00+. The point 
where the kinetic term changes its sign corresponds to 
the inflexion point of the scale factor (see Fig. 1). The 
branch At > represents a bouncing universe with a 
non-vanishing absolute minimum. 



B. Quintom 

We consider the standard quintom model with in- 
ternal metric Qik = diag(l,— 1) and assume that the 
scalar fields and (/)2 are coupled through the poten- 
tial V{(f>i, (j)2). So, the Fricdmann equation reads 



(43) 



where 01 represents the quintessence and 02 pl^'ys the 
role of a phantom field. The system is closed with the 
Klein-Gordon equations, one for each field. 

0i+3i7(/;i + 90,y(0i,02) =0, (44) 
02 +3if 02- 9^.^01 >2 ) = 0. (45) 

We propose a potential which can be written as follows. 



V^(01,02)=yo + |0? + ?02, 



(46) 



where a and (3 are constant and Vb > 0. It can mimics 
a simple model containing a mixture of dark matter and 
dark energy where the former is the clustering component 
and the latter give rise to the current acceleration. With 
the help of the Eq. PS|) we find that 



(47) 



where we have assumed separability of the potential 
F(0i,02) = V(0i) + V^(02). Thus, Eqs. (glD-dini) are 
easily integrated 



01 = 001 h = 



•2 = <P02 a 



-3/(1-/3) 



(48) 



whit 001 and 0o2 integration constants. Combining Eqs. 
(|43)) . (|46)) . (|48|) . the Fricdmann equation becomes 

~^°+2 a6/(i+-) +2 a6/(i-/3) ' ^^^> 

On the other hand, the effective barotropic index 7 = 
{p + p)/p of this mixture 

12 -6/{l+a) _ 12 -6/(1-/3) 

^ , (50) 



shows the occurrence of a smooth transition from the 
no phantom regime to the phantom one. This tran- 
sition is achieved when j{aph) = 0, that is at aph = 

(0O2/0Ol)(^-'^»^+'^'/^(^+«. 

In the case a = 1 and /? = 3, we associate the second 
term of the Fricdmann equation (|49|) with an attractive 
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component pm = 4''oi/o?i which includes both baryonic 
and nonbaryonic matter, having the equation of state 
p„j w 0, while the remaining two terms constitute the 



dark energy component pde = Vb 



The Eq. 



with 002 = ^ ^^^d (f^i = Vb/4 is solved by changing the 
variable to a 
± 



v^/'^. It turns into the equation Av 
2w^). whose solution gives the scale factor 



^/2 



tan ■ 



2/3 



(51) 



where At = t — to. Besides, the scalar fields (pi and 
are obtained after integrating Eq. 



V3 



cos 




(52) 
(53) 



The solution (j51|) represents a universe starting as it were 
dust dominated by the scalar field 0i with a « At^^^. 
After that it expands and ends in a big rip singularity at 
finite time. The big rip is reached where 02 blows up. 

Finally using the Eq. ([48|) we are able to obtain the 
potential (|46l) 



Vn r 



3e 



(54) 



as a sum of two separate potentials, each one depending 
on separate fields V^(0i, '?!>2)=^i(0i) + ^2(02)- This is in 
agreement with the precursor quintoni model proposed 
in the literature [2^ . 

Now, we investigate the degenerate case P = —a with 
CT > 0. It gives rise to singular and bouncing solutions 
with a final de Sitter behavior. The Friedmann equation 
(l49l) reads 



3i7^ = Va 



1 



^01 



'/'02 



a 1+° 



(55) 



Here there is a degeneration because the contribution 
of both fields arc proportional. When (p^-y > 0q2, we 
recover the ACDM cosmological model for cr = 1 and 
cosmic string for <t = 2. When 0§]^ < 0§2j the solution 
bounces where the total energy density vanishes. Solving 
the Friedmann equation ([55|) , we obtain the solutions 



\fh sinh ujtd 



-6 cosh a; At 



6 > 0, 
h < 0, 



(56) 
(57) 



where = 3Vo/{l + a)^ and b = (ct + l)((/)2^ - 022)/2Vb. 
The scale factor ((56|) evolves like a (x t^'^+^y^ near the 
singularity, having an inflationary phase for cr > 2 and 
ending in a de Sitter stage. The bouncing solution ([FT)) 
begins and ends with de Sitter phases. From Eqs. 



((56)) - ((57|) . after integrating, we find the fields in terms of 
the cosmic time 



yoi<? 



(58) 
(59) 



with 



1 , , wAi 
Intanh — - — , 



:tan-ie"^*. 



(60) 

Hence, coming back to the Eq. (|46|) we get the potential 
corresponding to cither case, quintessence or phantom 
dominated models 



V^Vo 
V^Vo 



1 H sinh^ Vbbjd)'^ 

1 + cr 



1 



1 + 0- 



-bLU(j) 



(61) 
(62) 



In this models there is no transition from non-phantom 
to phantom behavior because the barotropic index has a 
definite sign. 



C. 3-dimensional internal space 

Now, we study the case of a multi-scalar field model 
containing three fields with the following internal metric 



Q 



ik 




(63) 



It represents a scalar field configuration where the dimen- 
sion of the internal space associated with the quintessence 
sector Nq exceeds the dimension of the respective phan- 
tom sector Npi^ , so Nq > Np^ ■ The case where the phan- 
tom dominates over quintessence components Np^ > Ng 
can be obtained from the latter by changing the sign of 
the metric Qik — > ~Qik- 

We focus on the case where the multi-scalar field is 
driven by the potential 

1/(01,02,03) = Vb + ^0^-^01, (64) 



so the Friedmann equation reads 



3ff2 = Va + 



§0 

n3 



n3 ■ 



(65) 



The four terms of the total energy density can be ar- 
ranged as a two-fluid mixture, with positive energy den- 
sities, that are conserved separately 



Pb = Pde 

Pb + iHpb = 0, 



Pde 

20?o - -^ioa^ 



^ _ ^ (661 



'AH^dePde 



0, (67) 
(68) 
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and < 4Vb(/'fo- have identified pf, and pde with 
the baryonic and dark energy densities. The latter com- 
ponent comprises vacuum energy density Vq, a stiff fluid 
term ps = (fio/a^ and a kind of pressurless perfect fluid 
with negative energy density po = — 03o/a^- The last 
term mimics the negative part of the classical Dirac Field 
[2?! . Hence, the internal space permits to incorporate the 
negative part of the classical Dirac Field as a source of 
Einstein equation in a natural way. 

The general solution of the Friedmann equation (p5|) 
takes the following form 




cosh " 



1 



(69) 



where the initial singularity was fixed at t = 0. The dark 
energy crosses the phantom divide at Cc = {24>\q/ (jy^QY^^ 
where it reaches its minimum value pdec ~ ^'zol ^Vio- 



w(a) 




Figure 2: The state parameter Wde as function of the scale 
factor. The dashed line includes the Dirac term and the solid 
line does not include it. 

Basically, the dark energy crosses the phantom divide 
due to the term 0§o [23|. Concerning the nature of this 
parameter, note that the dominant energy condition is 
violated in this cosmological scenario thanks to the term 
(/)|q. When it is absent the dark energy state parame- 
ter Wde always remains in the range [—1, 1] (sec Fig. 3), 
implying that the dark energy component satisfies the 
dominant energy condition. 



V. SUMMARY AND CONCLUSION 

We have extended the quintom model by introduc- 
ing a multi-scalar field configuration along with a A^- 
dimcnsional internal space endowed with an constant in- 



ternal metric and given the lagrangian formulation of this 
new model. This formulation enriches the duality be- 
tween contracting and expanding cosmologies. We have 
shown that power law solutions arc asymptotically stable 
when the multi-scalar field is driven by the exponential 
potential with the same internal symmetry than the ki- 
netic energy term, generalizing similar results obtained 
for the quintessence cosmology. 

We have found the general solution of the nonlinear 
Einstein-Klein-Gordon equations for the exponential po- 
tential V = Voe~*^'^, extending previous works with 
standard scalar fields [l^ and k-essence [1^. The rich- 
ness introduced by the A'^-dimensional internal space and 
the internal metric Qik leads, in a natural way, to new 
singular solutions with a future big rip. For instance, in 
the explicit solution case, the scale factor blows up at 
finite time for < 0, exhibiting a superaccelerated be- 
havior and ending in a final big rip. It has a resemblance 
with the big rip behavior in the phantom universe re- 
ported recently by Wei [2^ for the potential V ~ Vb$^^. 

We have shown that the quintom model with a sepa- 
rable potential can be interpreted as a system of coupled 
fiuids. This framework contains the ACDM model, con- 
duces to bouncing universes and mimics the cosmic string 
cosmological model with cosmological constant. 

To gain insight into the internal space structure we 
have studied a model with a metric Qik, representing a 
configuration where the dimension of the internal space 
associated with the quintessence sector Nq exceeds the 
dimension of the respective phantom sector Npf^. Choos- 
ing (A^g, Nph) = (2, 1) we have accommodated a 3-scalar 
field with the following three components, stiff fiuid, dust 
and a dust with negative energy density. Then, the en- 
ergy density of the dark energy, composed with the lat- 
ter component, stiff fluid and vacuum energy crosses the 
phantom divide. Comparing with the Ref. [27| we see that 
the dust component with negative energy density plays 
the same role that the negative part of a classical Dirac 
Field. 
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The constant fca can be set equal to the unity scaling the 
time variable. Also, we choose k2 = — tq, then v = At = 



T - To. 

Without loss of generality we can choose 0oi 
where (f)o is an arbitrary constant. 



ilog\(t>o\ 



